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Abstract 

We study the Hamiltonian formulation of Plebanski theory in both the Euclidean 
and Lorentzian cases. A careful analysis of the constraints shows that the system 
is non regular, i.e. the rank of the Dirac matrix is non-constant on the non-reduced 
phase space. We identify the gravitational and topological sectors which are regu- 
lar sub-spaces of the non-reduced phase space. The theory can be restricted to the 
regular subspace which contains the gravitational sector. We explicitly identify 
first and second class constraints in this case. We compute the determinant of 
the Dirac matrix and the natural measure for the path integral of the Plebanski 
theory (restricted to the gravitational sector). This measure is the analogue of 
the Leutwyler-Fradkin-Vilkovisky measure of quantum gravity. 



I. Introduction 

Plebanski theory [T| is a 4-dimensional BF theory with an additional field which 
forces the B field to satisfy the simplicity constraint. By Plebanski theory we 
mean here the real formulation described by the real action defined below (|13|) 
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and not the complex formulation whose properties have been investigated in [21 El- 
It contains, as a particular sector, 4-dimensional pure gravity and is therefore an 
interesting field theory. The quantum properties of this field theory are, however, 
largely unknown. One important line of study aims at discretizing this quantum 
field theory with the tools of lattice gauge field theory leading to spin-foam models. 
Although spin-foam models have been the subject of numerous works over these 
last years (see the introduction jlj and the review 0), central issues are not 
understood and the technical tools needed to address these central questions need 
still to be developed. In particular we had in mind two pressing questions when 
beginning this work: 

-is it possible to compute from first principle the weight of the faces, edges and 
vertices in the spin foam model description of Plebanski theory? 

-can we see the appearance of quantum groups in Plebanski theory with cos- 
mological constant? 

The following work is a study of the Hamiltonian description of Plebanski 
theory. In particular we want to address the following problems: 

-computation of the Liouville measure in the path integral expressed in term 
of the original variables of Plebanski theory. This could be a first step for under- 
standing how to fix the measure of spin foam models. 

-computation of the Dirac bracket of all the fields once all second class con- 
straints have been taken into account. In particular one could have expected to 
find some brackets of quadratic type which are the signature of quantum groups. 
Unfortunately the Dirac brackets of the variables are complicated and no hint of 
quantum groups have appeared yet using this method. 

Our paper is organized as follows. In the next section (section II), we briefly 
recall the Plebanski description of gravity as a constrained BF theory. We then 
give the Hamiltonian analysis of the theory (primary and secondary constraints, 
first and second class constraints) (section III). We finally derive the path integral 
measure (section IV). 

II. Gravity as a Constrained Topological Field 
Theory 

This section is devoted to briefly reviewing the classical Lagrangian analysis of 
BF and Plebanski theories, and to presenting the notations that we are using 
throughout the work. Let M. be a four dimensional oriented smooth manifold Ai 
and the space of one-forms defined on the manifold M.. A local coordinate 

system will be denoted by x^, with \i 6 {0, 1, 2, 3}. 
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II. 1. BF action with cosmological term 
II. 1.1. Dynamics 

For any Lie group G one can define a BF theory which is a topological quantum 
field theory [HI El IE] • in the sequel we will only be interested in the case where the 
group G = SO(rj) with r\ = diag(a 2 , 1, 1, 1) and a = 1 (resp. a = i). We will coin 
the name Riemannian (resp. Lorentzian) 5F-theory in these cases. As usual, 
we will denote g = so{rj) the Lie algebra of the group whose basic properties 
are recalled in the appendix A. By convention, given two g-valued one forms 
a = a^dx^ and (3 = ft^dx' 1 , their commutator is defined to be the g- valued two- 
form [a, (3] = [otp, j3 u \dx >1 Adx u whereas the exterior product of these forms is given 
by a A /3 = -[a, /?]. We will also make an extensive use of the symmetrization 
notation: = -(7^ + j vft ) for any matrix 7. 

The action defining .BF-theory with a cosmological constant A is a functional 
Sbf[A,B] of a g-valued two-form B = ^B^ u dx^ A dx u , the 5-field, and a g- 
connection A = A^dx* 1 given by: 

S BF [A, B) = ± jj A x e^ pa ^< B^, Fp a > +| ~< B^, B pa >-j (1) 

where F = dA + A A A = ^F^dx^ A dx v is the curvature of the connection A, 
e iivp<y j s ^ e Levi-Civita symbol and the two independent Killing forms on the Lie 
algebra, < -, • > and -<•,•>», are precisely defined in the appendix A (see (JH2J))- 
As explained in the next section, the choice of the second Killing form in the 
cosmological term is made so that the full Plebanski action, obtained by adding 
the so-called simplicity constraints to (JTJ), is equivalent (in the appropriate sector) 
to ordinary gravity. 

Although this theory is well-known to be a topological field theory (see [1] for 
example), it still contains many unsettled questions such as the evaluation of the 
expectation values of observables (see jH] as well as more recent articles on this 
subject.) 

The equations of motion of BF action with respect to the fields B and A are 
respectively given by: 

t(F) + A B = and DB = dB + [A, B] = . (2) 

In these equations, 1 denotes the Hodge map defined in the appendix A. The 
situation is therefore quite different in the cases A = and A / 0. 

The B- field is locally covariantly constant but in the case where A 7^ 0, is 
completely determined once the connection A is given. Furthermore, in that case, 
the second equation (0) is an immediate consequence of the first one and the 
Bianchi identity. Any solution of the equations of motion is therefore given by a 
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couple (A, B) with A any connection and B given by B = —A~ 1 i(F). When the 
cosmological constant is zero, any solution of the equations of motion is given by 
a couple (A, B) with A any flat connection and B covariantly constant. Moreover, 
as the connection A is flat, the 5-field can be locally written as B = Dui for some 
0- valued one- form uj. 

II. 1.2 Gauge symmetries 

We now discuss the symmetries of the FF-action. First of all, we have the group 
of gauge transformations Q = C°°(Ai, G) whose action on the fields is given by: 

V geG,A 9 = g- 1 Ag + g- 1 dg and B 9 = g~ l Bg . (3) 

The invariance of the action (0) under the group of gauge transformations is 
an immediate consequence of the invariance of the Killing forms ([82)1 . Besides 
these familiar gauge transformations, .BF-theory admits another set of symmetries 
T = VL l (M) (gi q whose action on the fields depends explicitly on the cosmological 
constant as follows: 

VcjgT, A w = A-2Al(u) andF w = B + Du - 2Al([uj,uj}) . (4) 

The proof of the FF-action invariance under these transformations is a straight- 
forward computation and uses, as a central tool, the property (J%3*J) and the Killing 
forms invariance. 

The symmetry group S of the action (JTJ is obtaining by combining (jSJ) and 
(j3J), i- e -, is the semidirect product S = Q >< T where the action on the fields is 
given by (A, B) < (g,co) = {{A 9 Y , (B 9 Y). The composition law between two 
elements {g\,U\) and (c/2,^2) of S is given by: 

(01,^1) (02,^2) = (#i#2,w 2 + g 2 1 Lu 1 g 2 ) ■ (5) 

At this point, it is interesting to note that, in contrast to the three-dimensional 
case P|, the symmetry group structure (JSJ) is the same for any value of A. The 
cosmological constant appears only in the expression of the action of the normal 
subgroup T on the fields (J2J). 

The infinitesimal expression of the symmetries easily follows from (JSJ UJ) and 
reads : 

5 {u , T )A = Du - 2Ai(r) and S^B = B + [u, B] + Dt , (6) 

where u G C°°(M,g) and r G SV-(M) <8> 0. 

Moreover, being given by the integral of a 4-form, the BF-action is also in- 
variant under the group of diffeomorphisms of Ai. In the case of Riemannian or 
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Lorentzian BF theory, there is an interplay between the action of diffeomorphisms 
and the action of the symmetry group S. To see that this is indeed the case, 
let us compute the Lie derivative of the fields A and B along a vector field 
£ = ^dp. An easy calculation shows that these Lie derivatives are related to the 
infinitesimal symmetries JB} as follows: 

where u = ^A^ and r = ^ p B P pdx^. As a result, if the fields A and B satisfy 
the equations of motion (j2J), their transformations under diffeomorphisms can be 
expressed only in terms of the symmetries ©, which form a complete set in the 
sense of [TU] . 

It follows from the form of the gauge transformations that all the solutions 
are locally gauge-related on-shell. When A ^ 0, two solutions of the equations of 
motion are related by a symmetry cu. When A = 0, two solutions of the equation 
of motion are such that the respective connections are flat, so they are locally 
pure gauge and therefore are related by a gauge transformation. The respective 
B fields are covariantly constant, therefore locally there exists uj which relates the 
two. 

II. 1.3. Hamiltonian analysis 

Let us end up this brief review on BF theory by recalling basic facts (HI El about 
the Hamiltonian analysis (see ^1] for a recent complete treatment). 

For this purpose, we assume that spacetime has the structure M. = K. x £ 
where the real line K. represents time and S is an oriented smooth 3-manifold 
representing space. We will use a chart adapted to this foliation: x° is a 
coordinate on K. whereas x l ,i = 1,2,3 are coordinates on S. As usual we denote 

e ijk _ e 0ijfc_ 

The action (Q) shows that the spatial components of the connection Af and 
e^ k B^ k b f] ab = Il( e)a are canonically conjugate fields, i.e. 

{Af )a (x) ,%, )b (y)} = 5 b a 5 e , £ ,5f5(x,y). (9) 

Moreover, a straightforward calculation shows that the canonical Hamiltonian 
H BF [A, B] reads: 

H BF [A, B] = - [ d 3 x (e l ' k < B 0l , F jk > + < A , > + A -< B w , IT y) . (10) 
Jt, 



(<V)A)<?" + §J (7) 

S S 

{8{u,t)B)%" - epvpat? x . {E)a , (8) 



SAP 
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The components Aq )a of the connection and the components a of the 5-field 
appear to be non-dynamical variables: they are Lagrange multipliers which impose 
the following primary constraints: 

e^(F ife )+Air = (11) 

day = o. (12) 

These constraints are obviously contained in the Lagrangian equations 0, are 
first class and generate the infinitesimal symmetries. Transformation laws of the 
dynamical variables (A'f' a , rP (e)a ) under these symmetries are given by © with 
parameters u and r such that tq = 0. 



II. 2. Plebanski action 

Plebanski action pQ is a constrained Riemannian or Lorentzian BF theory, in 
which the simplicity constraint on the B field written explicitly in below is 
dynamically enforced through the introduction of Lagrange multipliers. 

The Plebanski action is a functional Spi[A, B, ip\ of the 5-field, the connection 
A and Lagrange multipliers <p. Actually, there exists two different but "equiva- 
lent" 5 formulations for this action ^2] depending on the choice of the Lagrange 
multipliers type: ip can be a density tensor field with space-time indices only or 
a field with Lie-algebra indices (internal) only. The latter formulation is more 
customary in the literature on spin foams but we will use the former which turns 
out to be more convenient for the achievement of an Hamiltonian analysis. The 
SO(rj) Plebanski action is given by: 

S pl [A, B, ip] = Sbf [A, B] + - [ d A x -< B pa y . (13) 

4 JM 

In this expression, ip^*** = c^IHIH j s a density tensor, symmetric under the 
exchange of the pair [fiis] with [pa] and satisfying the tracelessness condition: 
^tiupa ( P >1Uf " 7 — 0. Therefore it contains 20 independent components. Varying the 
action (fTB^) with respect to tp imposes the following 20 algebraic equations on the 
36 components of the 5-field: 

V 

B^ u , Bp a r" — £fj,upij ( 1 '1 ) 

where V = e M ^ pcr -< B^ u , B pa >- . In the sequel, V will be called the four-dimensional 
volume. When V is not vanishing it is said that the non degeneracy condition is 
satisfied. The equations (|14|) are known as "simplicity constraints" and their 

5 To be more precise, the two formulations are shown to be equivalent on shell when the non 
degeneracy condition is satisfied. 
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properties and solutions have been studied and classified when the non degener- 
acy condition is satisfied ^3 El E2] • In particular, it was shown that they are 
equivalent to the following system of equations: 

e ^ B U B KL = V e j JKL (15) 

Moreover, simplicity conditions and the non degeneracy conditions are simultane- 
ously fullfilled if and only if there exists a non degenerate cotetrad field e 1 = e^dx^ 
such that the 5-field can be written in one of the following forms: 

topological sector (I±) : B IJ = ± e 1 A e J (16) 
gravitational sector (II±) : B IJ = i ^ U kl e K A e L . (17) 

We now come to the relation between Plebanski theory and Palatini formulation of 
general relativity. If one substitutes the solutions of the sector (Hi) (fT7|) into the 
Plebanski action, one obtains ±<S±^[yl, e], where S^fAe] is the Palatini action 
of general relativity with cosmological constant A : 

S?[A } e] = i Jj'x e IJKL ej e J u (f* l + ± ef ety . (18) 

This result gives an a posteriori justification of the choice for the Killing form 
defining the cosmological term of the £?F-action (JTJ). The topological sector how- 
ever does not have any physical interpretation. 

Note that when the non degeneracy condition is fulfilled the different set of 
solutions of the simplicity conditions are disjoint. To see that this is the case, the 
Urbantke metrics g$, e G {+1,-1}, defined as [To] : 
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o 

(e) — p , ^«/?7<5 p>^ a D( £ ) 6 D( £ ) c (-\Q\ 

fiv ~ gy t abc t D lict fi-y n ' 5u • \ Li) J 

are a useful tool. If the S-field satisfies the simplicity conditions, their expressions 
simplify drastically and are given in term of the cotetrad by the formulae: 

topological sector (Ii) : g^J — ± e rju e J v (20) 
gravitational sector (Hi:) ■ gff = ± rju e J v . (21) 

As a result, a non-degenerate .B-field solution of the simplicity constraints cannot 
belong to both topological and gravitational sectors. Moreover, given a 5-field 
in the gravitational or topological sector, the signature of its associated Urbantke 
metrics completely fixes its belonging to the (±) sectors. The four sectors are 
therefore disjoint. 
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III. Hamiltonian analysis of Plebanski theory 



We now study SO(rj) Plebanski theory in the Hamiltonian framework. For this 
purpose, we will assume that M. = R x E where E is a smooth, oriented three- 
manifold (without boundary) which will play the role of space. 

The first step consists in deriving the Hamiltonian and the constraints. One 
must then split the constraints into first class constraints and second class con- 
straints. A notable feature of the Plebanski action is that it does not define a 
regular Hamiltonian system in the sense that the rank of the Dirac matrix de- 
pends on which type of sectors the S-field belongs to. 

We provide a set of holonomic constraints which replace the simplicity con- 
straints and eliminate the topological sector. When the S-field is non holonomi- 
cally enforced to belong to the gravitational sector 11+ we can compute the Dirac 
matrix determinant in closed form. This will be used in section IV to formulate a 
path integral quantization of Plebanski theory in this sector. 

III.l. Constraint analysis 

The Plebanski action (jlHj) can be written, after an integration by parts, as: 

S Pl = f dt J d 3 x(e mnp < B mn , d A p > +e mn ? < B 0p , F mn > + 
< A , D p B mn > e mn ? + A -< B 0p , B mn y e mn ? + 

(p 0m0n -< B 0m , B 0n y +^ 0mri P -<; B 0m , B np y +-(p mnrs -< B mn , B rs y). (22) 

From this expression we see that the fields Aq, Bo m , tp are non dynamical since their 
time derivatives do not occur in the Lagrangian. However the fact that the i?om's 
appear quadratically prevents us from treating them as Lagrange multipliers. To 
handle this problem, we decided to add to the Lagrangian the term 

<P i ,9 0J B 0l >-< /Ul ,P i > 

where P l and fii are new independent variables. This does not change the dynam- 
ics since the new fields are auxiliary fields, the elimination of which yields back 
the original action (by varying with respect to ^ and P l , one gets P l = and 
d()B i — ^ = 0). This modified action with the new fields is our starting point. 

III. 1.1. Conjugate momenta and primary constraints 

With the new variables, the Plebanski action is in first order form: 

S$[Ai, Hi, P h B 0i , A , m, if] = J dt Jd 3 x (< n\ doA > + <P\ d B 0i > -#)(23) 
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where we have made the change of variables 



IP = e mm B mn (24) 
and where the Hamiltonian H reads: 

-H = e mnp < B 0pj F mn > + <A , D P U P > +A -< B 0p , IP y + 

_|_^0m0n _^ B Qm ^ £> Qn y _|__^,0 m "P B Qm , IP >~ € npr + 



+^;<p m " r 'e mni e r8j ■< IP, ffV-< Mi , P l > . (25) 



We will denote V =-< .BojjIP y and assume that the non degeneracy condition 
holds i.e V 7^ 0. We conclude immediately that the only dynamical fields are the 
spatial components of the connection A4, their conjugate momenta IP as well as the 
Boi and their momenta P l . The other fields Aq, fa and <p do not appear with time 
derivatives. Moreover, the action (}2*3*j) is linear in them. They can, therefore, be 
considered as Lagrange multipliers. Thus, we do not associate conjugate momenta 
to these variables in order not to complicate unecessarily the canonical analysis 
and we start with the phase space of the A^s, the IP's, the B^s and the P*'s with 
following non-zero Poisson brackets : 

{A?"(x),%, )b (y)} = 8 a b 8 £>£ , 5} 6(x,y) , (26) 
W (x) , P? e , )b (y)} = 5 a b 5 £ , £ ,5]5(x,y). (27) 

The "primary constraints" are given by 

P[e)a « (28) 

T = D P W « (29) 
^B,B) t] =^B^B 0j y^ 0, (30) 

1 
3 

$(n, U) ij =-< IP, IF ' y « . (32) 



$(5, n)j =-! So,, n'y -77 5j v « , (3i) 



The constraint ()28|) is obtained by varying with respect to fa, constraint (}2"9"|) is 
obtained by varying with respect to A and the simplicity constraints (|30J31J32J) 
are obtained by varying with respect to ip. The time evolution of any phase space 
function / is given by / = {/, H}. 

III. 1.1. Secondary constraints 

We now study the secondary constraints arising from the conservation of the 
primary ones. 
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Conservation of the constraints (|28|) under time evolution implies the following 
conditions: 

pi = _ e Hk F]k + aA + 2aip mj l(Bqj) + a^ k e ]kl t{U l ) w . (33) 

These 18 equations involve the 6 components c/j * - 7 and the 8 components tp 01 ^. 
Due to the invertibility of V and the expressions of the primary simplicity con- 
straints, 14 of these equations completely fix the value of <^ 0l0 - ? , Lp°^ k whereas the 
4 remaining ones are secondary constraints given by: 

C = < P\ B 0l > = e^ k < B 0l , F jk > +A -< IT, B 0i y « , (34) 
d = e ijk < P\ n fc > = < F vv IF > « . (35) 

Conservation of the constraint (|29|) does not generate any new secondary con- 
straint. To see this it is useful (and permissible) to redefine the constraint T as 
follows: 

T = T-[P i ,B oi ]^0. (36) 

The constraint T generates infinitesimal gauge transformations of type Q on all 
the dynamical variables. As a result, the constraints T are conserved under time 
evolution because H is invariant under these transformations. 

Before studying the conservation of the simplicity constraints we prove some 
relations implied by the non degeneracy condition. We denote by K, L the 6x6 
matrices defined by 

K {e)a _ R (e)o K (e)a _ ™ j 3+j _ 3e ( £ ) a ■ 3e { 

A i — a 0i i A 3+j - LL (e)ai ■ L ( £ )a ~ Oi > L (e)a ~ ii (e)a 1 ~ L , •••) 6 - K 6 ' ) 

The non degeneracy condition and the simplicity constraints imply that L is the 
weak left inverse of K. This implies that L is also the right weak inverse of K 
which reads: 

E n «. 5 S?' + K'Xi* - ^<W- (38) 

k 

First, we consider the set of constraints ()30I31|) whose time derivative is given 

by: 

^{B,B)ij = 2^v (i ,B 0j) yn 0, (39) 
$(£?, n)5 =^ /a, IF y +2e imn -< B 0j ,D m B 0n y + < A , [B 0j , IT] > 

\sU^ fi ki Tl k y +2e kmn -< B ok , D m B 0n y + < A , [B ok , U k ] >) « 0. (40) 



These 14 equations fix 14 independant linear combinations among the 18 multipli- 
ers fi\ £ ^ a and imply no further constraint on the dynamical variables. Indeed, let 
us introduce the variables v™ with //^ = Ylm=i v™Km ■ The set of constraints 
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flgmnD fix the combinations z/| +i + and v) - \5) Yf k=1 v\ where i, j £ {1, 2, 3}. 
From this analysis, we can reorganize the 18 constraints P l ~ in two subsets of 
constraints: 

k =< P\B 0i >^0 (41) 
Ki = e ijk < Pi,Tl k >^ (42) 

$(£?, P)» =< P\ B 0j > ~8) < P k , B ok >^ (43) 

$(P, Tl) ij =< P ( \ Tl j) >^ 0. (44) 

We will show later that kq , k% are the 4 first class constraints corresponding to the 
previous four Lagrange multipliers among the /i's that remain arbitrary. 

The study of the conservation of the remaining simplicity constraints ()32|) 
introduces, by contrast, new constraints in the theory. Indeed, vanishing of the 
time derivatives of these constraints imposes the following conditions: 

-i$(n,n) ij ' = T ij = e mn{i -< D m U j \B 0n ,y « . (45) 
This closes the analysis of the secondary constraints. 

III. 1.3. Conservation of secondary constraints 

In order to complete the Dirac algorithm we have to study the conservation of 
these secondary constraints. A direct computation shows that time derivative of 
the constraints (|45j) can be expressed as: 

fij = _ 2€ uv(i e m j)p ^ Du D m B 0p , B 0v y +e m »* -< [A 0j B 0n ],D m W ] >- + 

+ a 2 Ae mn {t < ni)j [Bom}Bon] > +e mn(i ^ ^jp)^ y +^0^ < Tp), [B 0k , B 0l ] > 

+±o*<p™ r 'e mnk e rtl < itfVW [W\B Gv ] >^ . (46) 

As a result, conservation of the constraints T 4 - 7 is a system of equations which 
will eventually fix some of the yet unfixed Lagrange multipliers ip klmn . This point 
depends crucially on invertibility properties of the following tensor: 

Quvij = a 2 < n (t> e «)rp [W\B 0r ] > . (47) 

We now show that G um i can be rewritten in terms of the Urbantke metrics. It 
follows that its properties depend on the sector under consideration. 

Lemma 1 The tensor G um i satisfies the identity: 

a 2 

Quvij w _ V ( £ ki(u e v)lj + e kj(u e v)n^ g {+) + g (-)y ^ 
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This implies that: 

-in the topological sector (I) the matrix G identically vanishes, 
-in the non topological sector (II) the matrix G is invertible and its determinant 
reads: 

det(G) « ~V\det{g kl )f (49) 

Proof: Using the formula defining the Urbantke metrics, and the formula ex- 
pressing Bij in term of IT, we get: 

glf « ^e abc e m{ir B^Uf E)b W {E)c . (50) 

Using this result we can now compute the right handside of (jUJ), which is equal 

to G uvij . 

We have previously seen that the topological phase (I) is characterized by 
the condition gffi + gj^/ = 0. This implies that G is zero in this sector. In the 
non-topological sector we have G uv ^ ~ ^-V(e fc ^ w e^^ + e k ^ u e v ^ h )gki- From this 
expression one easily computes the determinant det(G): 

det(G) « cV\det{g kl )f (51) 

where the constant c is evaluated to be c = —\ (go to a frame where g k i — Ski). □ 
In the sequel, we will denote by V s the spatial volume defined by V, 2 = 
\det(gki)\, so we have det(G) « cV 6 Vg. 

Because the matrix G is invertible in the non topological sector, the equations 
T %: > = fix the value of the 6 remaining Lagrange multipliers (p Mmn and do not im- 
pose any new constraints. By contrast, in the topological sector, these equations 
do not fix these Lagrange multipliers anymore but give tertiary constraints. The 
Plebanski theory turns out to define a non-regular Hamiltonian system, i.e. its 
constraint analysis is not uniform on the non-reduced phase space. In the topo- 
logical sector, we could continue the Dirac algorithm and we would observe that 
the set of first class constraints is larger than in the non topological sector. The 
gauge symmetries generated by this set of first class constraints would absorb all 
the local dynamical degrees of freedom. This justifies the term "topological" for 
this sector. In the sequel, we will exclusively concentrate on the non topological 
sector. 

It now remains to understand the conservation in time of the constraints Cq, Cj. 
By anticipating the discussion concerning separation of the constraints between 
first class and second class, we will modify the constraints Co,Cj by adding to 
them linear combination of constraints in such a way that the resulting constraints, 
denoted C^C^, are first class. 
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This is possible: let us define 



C[ = C i +<P j ,D i B 0j >-<A i ,T>-d j <P j ,B 0i > (52) 
= < diAj, IF > -dj < A, IF > + < P j , diB Gj > -dj < P j , B 0i > . 

We have verified by direct computation that C[ are indeed first class constraints 
(the only non trivial computations are the weak Poisson commuting property with 
respect to T* J and Co )• Consequently, they are conserved under time evolution. 
As to the constraint C' , we define it as being the total Hamiltonian with the 
values of the Lagrange multipliers fixed by the Dirac analysis: this is precisely 
solving the condition that C is first class with C' Q — C a linear combination of 
constraints. 

III. 2 First and Second Class constraints 
III. 2.1. Separation of constraints 

Once all the constraints have been determined, one must split them into first class 
and second class constraints. To achieve this (in the non topological sector), we 
reorganize the set of all the constraints as follows: 



Ko =< P\ B 0i (53) 

m = e ijk < P 3 ,U k >^ , (54) 

$(£, P)j =< P\ B 0J > ~5y < P k , B ok >^ (55) 

$(P, U) ij =< P ( \ U j) >^ (56) 

T = D p W-[P p ,B 0p }^0 (57) 

$(£,£);,- =^ £ i, #o, ^~ (58) 

$(p, n)j =^ B 0j , nV-^;v«o (59) 

$(n, u) ij =^ it, if >-« o (60) 

rpj = e mn(i ^ Dm n J ), B 0n (61) 

C[ =< F lv IF > - < P\ DiB 0j > + < Ai, T (62) 

C « (63) 



By a direct computation it is easy to check that Kq, k,i, Y^ a , C' , C[ are first class 
constraints. The remaining constraints turn out to be second class. To show that 
this is indeed the case, we compute their Dirac matrix (i.e., the matrix of the 
Poisson brackets of the second class constraints) and prove its invertibility. 

It appears that all the Poisson brackets of the remaining constraints are weakly 
ultralocal except {T lJ (x),T kl (y)} which contains non ultralocal term of 8'(x,y) 
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type. This can be taken care of (leading to simpler brackets) by redefining T tj as 
follows: 

T>(x)=T>(x) + J dyX^(x,yMU,U) kl (y) 

where X^(x,y) = — \{T ij (x), T k ' v (y)}G^ kl (y), so that all the Poisson brackets 
are now weakly ultralocal (i.e., contain only 5(x,y) and none of its derivatives). 
Thus, if we denote by \ a the set of redefined second class contraints, we now have 

{Xa(x),xp(y)} = A af3 (x)5(x,y). 



We shall exhibit the coefficients A a p. 

The Dirac matrix has the following structure: 



A a ' a 


®(B, B) mn 


$(n,p) rs 


HB,ny v 


HB,P)i 


$(n,n) w 


rpkl 


®(B, B) m i n i 





Ars 
m'n' 














$(n,p) rV 


Ar's' 
mn 













j^r's'kl 













B u ', j 

v'l 





jvru'kl 
v' 







_ R frs 


-B j ', u 

i'v 








Qi ki 


$(n,n) pV 

















£tp'q'kl 


rfk'V 





-M rsk ' 1 ' 


jvruk'l' 




_Qpqk'V 






where the tensor A, B, K, M, N, Q can be easily written down. In order to compute 
the determinant, only the values of A and B are necessary (besides the already 
known G) and we shall thus give explicitly here only those. We have: 

Wn'M = I^IA (64) 
Bii(x) = ±V(S:,5f-±5{5S) (65) 

Proposition 1 The determinant of the Dirac matrix is: 

D 2 = det A w det(A) 2 det(B) 2 det(G) 2 = 2 8 3" 28 V 4 V S 8 . (66) 

Proof: A developpement of the determinant of the Dirac Matrix along lines and 
columns give immediately the first equality. The second equality is obtained by 
pluging the explicit expression of det (A), det (B), det (G) which are given by: 

det(A) = (^) 6 V 6 ,det(5) = (i) 8 V 8 . (67) 

□ 

Because this determinant does not vanish, the constraints Xa are indeed all 
second class. 
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III. 2. 2. Dirac brackets 



Due to the simple structure of the Dirac matrix, the expressions of the cofactors are 
easy to obtain explicitly and appear to be polynomials in the dynamical variables 
and their derivatives. This implies that we can compute exactly the inverse of the 
Dirac matrix and we will denote by {, }d the Dirac bracket. We will not give the 
explicit expression of the Dirac bracket of the dynamical variables in the present 
work because no interesting structure seems to show up. In particular, we have 
seen no immediate hint for a quantum group structure. 

III. 2. 3. Counting degrees of freedom 

The counting of the degrees of freedom is direct: there are 72 dynamical variables 
(Ai, IT, B 0i , P l ), 40 second class constraints and 14 first class constraints. There- 
fore, the reduced phase space is 4 dimensional, which is correctly the number of 
degrees of freedom of general relativity. 

The 14 first class constraints k , Kj, C , T^ a are related to the gauge sym- 
metries of the theory. As was pointed out above, the constraint functions Y( e ) a 
generate infinitesimal gauge transformations of type Q. As to the constraints 
kq, K{, Cq, C'i, they are related to the spacetime diffeomorphisms. The primary 
constraints Ko = 0, Ki = appear because we treat the components B 0i , which 
contain the lapse N and the shift N\ as dynamical variables (so, the constraints 
K Q = 0, Ki = are equivalent to the constraints 7rjv = 0, n N i = 0). The secondary 
constraints C' Q — 0, C[ — are equivalent to the Hamiltonian and momentum 
constraints. 

Note that as a consequence of the constraints, the Hamiltonian H weakly 
vanishes. This was somehow to be expected, since the theory is diffeomorphism- 
invariant. 

IV. Path Integral Measure 

This section aims at defining the path integral of the Plebanski theory restricted 
to the gravitational sector. Many difficulties have to be overcome: 

1. to formulate a formal path integral taking into account the second class 
constraints 

2. to apply a BRST procedure to gauge fix the first class constraints 

3. to regularize and/or to discretize these path integrals 

4. to "renormalize" these amplitudes 

Our (modest) contribution will cover only the first point. We will define the path 
integral measure in term of the Liouville measure on the phase space taking into 
account the second class constraints. We will then apply the general method 
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of JHI to rewrite this path integral measure as a path integral on the original 
configuration variables of the Lagrangian formalism. 

We first give a brief introduction to the path integral method for pure second- 
class constrained systems [T71 [TBI [TO] , which is rather direct. We start with a 
hamiltonian system with 2n coordinates (x 7 )/ on a symplectic phase space (V,uS) 
and a set of 2m regular and irreducible second class constraints (Xa)a ■ The 
reduced phase space (C,u>d) is the submanifold C = {p G V\xa(p) = 0,Va = 
1, ■ • ■ , 2m} endowed with the Dirac symplectic structure u>d- 

The path integral is formally defined as a sum over all paths satisfying the 
constraints and the integral measure is given by the Liouville measure associated to 
uj d . To be more precise, let us introduce a complete set y l of 2{n — m) independent 
coordinates on C and we obtain the following definition of the path integral: 

Z = f \Vy\t)} \{ ^detu D (y(t))eMiS[y(t)}) . (68) 

In the previous formula, S[y(t)] is the action of the path y(t) and detun is the 
determinant of the symplectic structure. 

For our purpose, it will be useful to rewrite the path integral in terms of the 
original variables (x 7 )/. To that end, we perform a coordinate transformation 
of the non-reduced phase space V from the coordinates (x 7 )/ to the coordinates 
(y\Xa)i,a- To go further, we use the fact that it is always possible to find co- 
ordinates (y l )i of the submanifold C such that {y\x&} = 0. The local symplec- 
tic structure of V is therefore given by the tensor product uj = uj£> © uj x where 
{^^f)a,f3 is the inverse of the Dirac matrix A Q/ 3 = {XocXp} an d defines the two 
form uj x = {(jj'^)dxa A dxp- 

As a consequence, the equation relating the volume elements in each coordi- 
nates system is obtained using the Liouville measure conservation under coordi- 
nates transformations: 



FJrf^ = ^= \\d X a n<iy 4 . (69) 

From this result, one obtains the following expression of the path integral in term 
of the original variables x 1 : 

Z = ( [Vx{t)\ Y[{^detu{x{t)) y/detA(x(t))) 
J t 

2m 

II e*p(iS[x(t))) . (70) 



2m 

A 

a=l 



Thus, we see that the key point in writing the path integral is the computation 
of the Dirac matrix and the evaluation of its determinant. This task was carried 



16 



out above for the Plebanski theory in subsection III. 2.1 and the determinant was 
denoted D 2 . 

In the case at hand, there are also first class constraints. These can be gauge 
fixed by using 14 gauge fixing functions £ a chosen such that D\ = (det{^ a , C/3}d) 2 
does not vanish on the phase space (we collectively denote the first class constraints 
by ipa). Once the gauge is fixed, there are only second class constraints left and 
one can apply the above formula. This yields (for the Plebanski theory restricted 
to the gravitational sector) 

Zpi = JlVA^VU^VBo^VP 1 } 

exp (i j dt J d 3 x (< U\d Ai > + < P\d B w >)^j 

x i($(n,p))i($(B,B))i($(B 1 n))i($(n,n))i(T) (71) 

Recall that the Hamiltonian vanishes when the constraints are taken into account. 
The gauge fixing conditions £ a depend on the canonical variables and can also 
involve explicitly the time ("time dependent canonical gauge conditions"). For 
the subsequent discussion, it will be simpler to take gauge conditions that do not 
involve the potential Ai. This is equivalent to taking gauge conditions that involve 
the tetrads only (and their derivatives), but not the connection. The boundary 
conditions to be imposed on the paths that are summed over in the path integral 
should be dealt with along the lines of [T§] . 
We may re-express (fTTj) as 

Zpi = JiVA^VU^VBo^VP 1 } 

exp (i J dt jd*x {< n\d Ai > + < P\doB 0i > -Hjj 

x£M xGM 

x <y($(n, p))6($(b, b))5($(b, n))5($(n, u))S(f) (72) 

since H, defined as 

-H = e mn v < B 0p , F mn > +A H B 0p , W y, (73) 

vanishes on the constraint surface. 

Our aim is to rewrite this expression as a path integral over the original co- 
variant configuration variables with the original Plebanski action. To proceed, 
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we first integrate over the momenta P % . This can be done by using the following 
change of variables: 



[VP 1 ] = J[V$(B, P)\ P)} [Vkq] [VKi] (74) 



where J is the Jacobian of the linear transformation P % \— > ($(5, P), $(11, P), fCo ; K 
easily shown to be equal to J — (-< B ok , U k y)~ 9 = V~ 9 . As a result we obtain 



I J ! 



Z Pl = I [VA^VU^VBo,] 

exp (iJ dt J d * x (< n*,flbA > -Hjj 



J] V^C^*(Co)*(C*)*(T)*(e.) 

rGM 

x y^V~ 9 5(^(B, B))5($(B, n))<5($(n, n))5(f ). (75) 



We now apply the method of 16J to eliminate the secondary second class 
constraints T l K We denote Ci = ~ i^(n, n)^ as well as = T 1 ^ . The constraints 
Ci are the secondary constraints associated with Ci i n the sense that we have 
{Ci , H} = Ci ■ Introducing multipliers Ajj and pij to enforce the constraints, we 
can rewrite (|73j) as: 



Z Pl = J [VAiWV^WVB^VX^Vp^] 

exp \i j * dt J d 3 x (< IT, doAi >-H- XijC? - PijC^jj 

II v / ^^^(^)5(a)5(T)5(e Q )v^v- 9 5($(P,P)) ( 5(<i>(p,n)). 

We consider a canonical transformation generated by the function F = — Cy- 
linder this transformation the kinetic term < IT 1 , d^Ai > is left invariant 6 , as well 
as the path integral measure and also \jCi , whereas H is turned into 

H' = exp({F, .}){&) = H- l>;j a - 2p ijPkl G ijkl 

and (2 i s turned into 



6 The canonical transformation depends explicitly on the parameter py, so the kinetic term 
transforms with the explicit time derivative of the generator, S < IP, Qo-Aj >^ PijCi ■ Since 
the constraint Ci remains invariant (and enforced through the integration over Ay), the kinetic 
term is indeed invariant on the constraint surface. Alternatively, its variation in the non-reduced 
phase space can be absorbed through a redefinition of Xij. 
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(to obtain the latter results, we have used {£?', C 2 M } = -4G^ and {G^' fcZ , F} = 0). 
Finally, since only Ai transforms under the canonical transformation generated by 
F, the product of the ^-distributions in the path integral is invariant (we assume 
that the gauge conditions do not involve AA. As a result, we obtain 



Z Pl = J [DAiliVU^VBo^Vp^} 

exp (ijdtj d3x (< R^doAi > —H — 2 PljPkl G tjkl ^j 
J] >/5^^(CS)*(CjWT)^ 

x&M 

The integration over the variables pjj is gaussian and gives 

Z Pl = J [VAi][VW][VB Qi ]ex. V (i J dt Jd 3 x (< ir,$,A > -H 

II v ^i^<5(^)5(q)5(T)^a)^^^5($(s ! s))<5($(s,n))5($(n,n)). 

, e t/ \/det(G) 

We have thus succeeded in eliminating the constraints T^. 

To make contact with the Plebanski action, we now reintroduce the original 
configuration variables. Because II 1 = e^ k Bjk, the Jacobian of the transformation 
from IP to Bjk is a constant, therefore [DIP] = [DBij]. We can furthermore elimi- 
nate the Dirac distributions 5(T), 6(<&(B, B)),5($(B, II)), 5($(II, II)) by restoring 
the integration over A$ and ^ vpa '. Taking into account the known values of yfD^ 
and a/ det(G), we obtain the final expression: 

Z Pl = J [VA^}[VB^][V^}V 8 V S 2 ^P (iS Pl [A,B,tp}) 

J] y/D^SiC'MC'^). (76) 

x€M 

This is the formal expression of the gauge fixed measure of Plebanski theory in 
the non topological sector. One can be puzzled by the fact that the measure has 
an expression which is non covariant because it contains V^. This phenomenon is 
typical of derivative couplings and is completely analogous to the non covariant 
measure computed 40 years ago by H.Leutwyler j20] which reads: 

Vg = l[V- 5 V s 2 U d ^ ( 77 ) 

X /i<I^ 

A derivation of this measure from the Liouville measure was done by E.S.Fradkin 
and G.A.Vilkovisky j2TJ [22j as well careful analysis of the covariance of the 
correlation functions. See also [23] for the different comparisons of functional 
measure for quantum gravity. 
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V. Conclusion 



In the present work we have derived the hamiltonian analysis of the Plebanski 
theory. A notable feature of this theory is that it defines a non-regular theory, 
in the sense that the matrix of the Poisson brackets of the constraints does not 
have constant rank on the constraint surface. Non-regular theories used to be 
thought of as being somewhat pathological and of no interest to physical systems. 
In the last years, however, models have been constructed that lead to non-regular 
constrained Hamiltonian systems. Another example, besides Plebanski theory, is 
given by the higher- dimensional pure Chern-Simons theories (211123 • A recent 
analysis of non- regular systems may be found in [2B|. 

We have explicitly carried out the derivation of the Hamiltonian formalism 
in the non topological sector. In particular, we have obtained all the constraints 
and been able to disentangle them into first class and second class. We have also 
computed exactly the local measure, which has enabled us to give a formal path 
integral of Plebanski theory (in the non topological sector) in term of the original 
variables. This constitutes a first step in the direction of the understanding of 
Plebanski theory. 

Our work should be further extended in various directions. First, the above 
path integral was written in a canonical gauge. It would be useful to derive it in 
propagating gauges by BRST methods. 

Second, it would also be of interest to try to discretize the Plebanski theory 
using spin foam models. A central issue is the computation of the weights of the 
faces and edges which are still a matter of debate j23 I2H1 I2H] • The computation 
of the functional measure that we provide will hopefully give new insights in this 
intricate problem. 

Finally, one of our motivations was to try to uncover a quantum group struc- 
ture underlying the Plebanski theory. In euclidean BF theory with cosmological 
constant it is well known that quantum groups are hidden in this topological field 
theory [HDj. Up to now there is no clear picture of appearance of quantum groups 
(more precisely of Poisson Lie groups or Poisson Lie groupoids) as symmetries 
of Poisson brackets of the fields in this theory. Unfortunately, the form of the 
Dirac brackets that we can compute in Plebanski theory does not appear to be 
connected with such structures in any obvious way. The relation, if any, between 
quantum groups and Plebanski theory still remains to us a mystery as well as the 
relation between the quantum deformation of the Barrett-Crane model (31] and a 
discretization of Lorentzian Plebanski theory with cosmological constant. 
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Appendix 

Algebra so(rj): definition, notations and properties 

g = so(r)) is the real Lie algebra of the isometry group G = SO(rj) of the quadratic 
form rj = diag(cr 2 , 1, 1, 1) where a G {1,?}- We will use capital latin letters for 
internal vector indices I, J, - ■ ■ G {0, 1, 2, 3} and we will define the antisymmetic 
tensor e IJKL such that e 0123 = 1; the indices are lowered and raised with the metric 
V- 

One of the basis of so(rj), composed of rotation generators (J a )a=i,2,3 and boost 
generators (i^ a )a=i,2,3 has the following commutation relations: 

[J a, M = tab-h , [K a , K b ] = a 2 e ab c J c , [K a , J b ] = e ab c K c . (78) 

We have used the notation e abc = e° a bc- 

The real Lie algebra so(rj) is also a real form of so(4, C) = so(3, C) © so(3, C). 

It is convenient to introduce the basis (T(±) a ) a =i,2,3 defined from the previous 
one by T(±) a = \{J a =t o'Ka) whose generators realize two commuting copies of 
so(3), i.e.: 

[T(±) a ,T(±)b] = e ab T(±) C , [T(+)a, T(-)b] = . (79) 

We denote by p + and p~ the projectors associated to this direct sum. To any 
element f of so(4, C), we denote ^ = p ± (0 = ^ ±)a T (±)a . 

The star structure on so(4, C) selecting the real form reads: 

*(T(±)a) = — T (CT 2 ±)a . (80) 

As usual it is customary to define the Hodge linear map t = -{p + — p~) which 
satisfies i : so(rj) — > so(r]) and id, as well as: 

V^e so(rj) , [t(0,x] = (81) 

The space of real symmetric invariant bilinear forms on the Lie algebra so{rj) is 
a two-dimensional vector space; we choose two independent nondegenerate bilinear 
forms <■, ■> and -<•, •>- defined by: 

< T(e)a , T( e i) b >= 5 ab 5 ££ i , -< T(e)a,T( e >) b >-— OS 5 a b b~se' ■ (82) 

From the definition of t, we have: 

ve,x e so(rj) , -< l(£),x y = -< Z,l(x) >- = < Z,X > • (83) 

The vector representation LT of the Lie algebra so(rj) is given by: 

n(T (±)a ) /J = T& a = i(-^±a(7/V 7 -^ V J )) • (84) 



22 



To prove that the formulae ()84j) define a representation of the Lie algebra, it 
is sufficient to verify the following useful relations: 

Ttf aVj jT$ b = 5 e£ , Qe afc c Til ~ \v<*S IK ) ■ (85) 

Furthermore, the values of the quadratic Casimir tensors of so(r]) in the vector 
representation can be computed using the relation 

T H T ^ = \^ UKL + \(8 IK 6 JL - 6 IL 5 JK ) . (86) 

a 

From these quadratic formulae, it is easy to compute the components of any 
element £ in the Lie algebra in term of £ /J and one obtains: 

£(±)« = _ v ab^IJ ^ T KL (g7) 
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